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Abstract 

A study of Bose-Einstein correlations in pairs of identically charged pions produced in e + e~ 
annihilations at the Z° peak has been performed for the first time assuming a non-static emitting 
source. The results are based on the high statistics data obtained with the OPAL detector at 
LEP. The correlation functions have been analyzed in intervals of the average pair transverse 
momentum and of the pair rapidity, in order to study possible correlations between the pion 
production points and their momenta (position- momentum correlations). The Yano-Koonin 
and the Bertsch-Pratt parameterizations have been fitted to the measured correlation functions 
to estimate the geometrical parameters of the source as well as the velocity of the source 
elements with respect to the overall centre-of-mass frame. The source rapidity is found to 
scale approximately with the pair rapidity, and both the longitudinal and transverse source 
dimensions are found to decrease for increasing average pair transverse momenta. 
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1 Introduction 



The space-time evolution of a source emitting particles can be probed using intensity inter- 
ferometry. Bose-Einstein correlations (BECs) in pairs of identical bosons have been studied 
at different centre-of-mass energies and for different initial states (e + e~ [1], pp and pp [2], 
lepton-hadron [3], nucleus- nucleus collisions [1]). BECs manifest themselves as enhancements 
in the production of identical bosons which are close to one another in phase space. They can 
be analysed in terms of the correlation function 



where p\ and p 2 are the 4-momenta of the two bosons, p(p\,P2) is the density of the two 
identical bosons and Po(pi,P2) is the two-particle density in the absence of BECs (reference 
sample). From the experimental correlation function one can extract the dimension of the 
source element (frequently called correlation length or radius of the emitting source), i.e. the 
length of the region of homogeneity from which pions are emitted that have momenta similar 
enough to interfere and contribute to the correlation function. 

At LEP Bose-Einstein correlations were analysed extensively in Z° hadronic events [5-11]. 
Two-pion correlations were studied as a function of the relative 4- momentum q = [pi — P2) 
of the pair: C(p\,p2) = C(q). It was found that the radius of the emitting region, supposed 
spherical, is of the order of 1 fm and increases with the number of jets in the event [5]. No 
significant differences were observed in the source dimensions between the 7r 7T ± and the 7r°7r° 
systems [Bj; on the other hand, smaller radii were measured in K ± K ± and K°K° (K°K°) pairs 
compared with pion pairs [7j. Genuine three-pion BECs were also observed [8]. Up to fifth-order 
genuine correlations of identically charged pions were obtained by OPAL [9] , where BECs were 
shown to be an essential ingredient of the correlation scaling observed there, also for all-charged 
higher-order correlations. The hypothesis that the source is spherical was tested studying the 
correlations in terms of components of q: two- and three-dimensional analyses have shown that 
the pion emission region is elongated rather than spherical, with the longitudinal dimension, 
along the event thrust axis, larger than the transverse one [TUl [EE]. BECs were also studied in 
e + e~ — > W + W~ events: no evidence of correlations between pions originating from different 
W bosons was found [12] . 

All the results listed above were obtained under the hypothesis that the momentum distribu- 
tion of the emitted particles is homogeneous throughout the source elements, as would happen 
if the source is static. In the case of a dynamic, i.e. expanding, source, the dimension of the re- 
gions of homogeneity varies with the momentum of the emitted particles. The expansion leads 
to correlations between the space-time emission points and the particle 4-momenta (position- 
momentum correlations) which generate a dependence of the BEC radii on the pair momenta. 
In this case, the correlation function is expected to depend on the average 4-momentum of 
the pair K = (pi +p 2 )/2 in addition to the relative 4-momentum q: C(pi,p 2 ) = C(q,K) [13J, 
so that the measured radii correspond to regions of homogeneity in K, i.e. effective source 
elements of pairs with momentum K. 

Published investigations of the source dynamics in e + e~ collisions are available at energies 
lower than LEP's [TJ. A dependence of the source radii on different components of the 4- vector 
K has been observed in more complex systems such as the emission region created after a high- 
energy collision between heavy nuclei. In particular, source radii have been found to decrease for 



increasing pair transverse momenta k t (or, equivalently, transverse masses m% = Jk% + m£) [I]. 



C(pi,p 2 ) 




(1) 




4 



Hydrodynamical models for heavy ion collisions [H] explain this correlation in terms of an 
expansion of the source, due to collective flows generated by pressure gradients. A similar 
dependence of the size parameters on m t was measured in pp collisions [15J. Longitudinal 
position-momentum correlations can be expected in e + e~ annihilations as a consequence of 
string fragmentation [16] . Models based on different assumptions (the Heisenberg uncertainty 
principle [17], the generalized Bjorken- Gottfried hypothesis [18J ) predict radii decreasing with 
the transverse mass also for sources created in e + e~ collisions. 

In this paper, which continues a series of OPAL studies on BECs [TTJ] , a measurement of 
three-dimensional Bose-Einstein correlation functions is presented and the correlation functions 
are analyzed in order to measure their dependence on K and investigate potential dynamical 
features of the pion-emitting source created after an e + e~ annihilation at a centre-of-mass 
energy of about 91 GeV. 

2 Experimental procedure 

A detailed description of the OPAL detector can be found in [191 I2Q] . In the present analysis, 
we have used the same data sample, about 4.3 million multihadronic events from Z° decays, and 
have applied the following selection cuts on tracks and events, identical to the ones described 
in [10]. First, the event thrust axis was computed, using tracks with a minimum of 20 hits in 
the jet chamber, a minimum transverse momentum of 150 MeV and a maximum momentum 
of 65 GeV. Clusters in the electromagnetic calorimeter are used if their energies exceed 100 
MeV in the barrel or 200 MeV in the endcaps. Only events well contained in the detector were 
accepted, requiring Icos^™^! < 0.9, where thrust is the polar angle of the thrust axis with 
respect to the beam axis q Then, a set of cuts, specific to BEC analyses, were applied. Tracks 
were required to have a maximum momentum of 40 GeV and to originate from the interaction 
vertex. Electron-positron pairs from photon conversions were rejected. Events were selected 
if they contained a minimum number of five tracks and if they were reasonably balanced in 
charge, i.e. requiring \n^ h — n~ h \/(n^ h + n~ h ) < 0.4, where n^ h and n~ h are the number of positive 
and negative charge tracks, respectively. About 3.7 million events were left after all quality 
cuts. All charged particle tracks that passed the selections were used, the pion purity being 
approximately 90%. No corrections were applied for final state Coulomb interactions. All data 
and Monte Carlo distributions presented here are given at the detector level, i.e. not corrected 
for effects of detector acceptance and resolution. 

The correlations were measured as functions of two different sets of variables, components 
of the pair 4-momentum difference q in two different frames. 

The first set, (Qi, Qt sidc , <5t out ), was evaluated in the Longitudinally CoMoving System 
(LCMS) [21]. For each pion pair, the LCMS is the frame, moving along the thrust axis, in 
which the sum of the two particle momenta, p={pi+P2)i lies in the plane perpendicular to the 
event thrust axis. The momentum difference of the pair, Q = {pi — P2) is resolved into the 
moduli of the transverse component, Q t , and of the longitudinal component, Qe, where the 
longitudinal (£) direction coincides with the thrust axis. Qt may in turn be resolved into "out" , 
Qtout) ancl "side", Qt sidc , components 

Qt = Qt out d + Q tsii J (2) 

x The coordinate system is denned so that z is the coordinate parallel to the e + and e _ beams, with positive 
direction along the e _ beam; r is the coordinate normal to the beam axis, is the azimuthal angle and 9 is the 
polar angle with respect to +z. 
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where 6 and s are unit vectors in the plane perpendicular to the thrust direction, such that 
p = pd defines the "out" direction and s = I x 6 defines the "side" direction. It can be shown [22J 
that, in the LCMS, the components Qt side and Qg reflect only the difference in emission space 
of the two pions, while Qt out depends on the difference in emission time as well. 

The second set, (g t ,g£,g ), was evaluated in the event centre-of-mass (CMS) frame. For 
each event, two hemispheres are defined by the plane perpendicular to the thrust axis. Each 
pair is then associated to the hemisphere containing the vector sum of the three-momenta. 
The pair 4- momentum difference q is resolved into the energy difference go — (Ei ~ E 2 ) and 
the 3-momentum difference q = —p^). The vector g is further decomposed into q t and qg, 
the transverse and longitudinal components, respectively, with respect to the thrust axis. In 
each pair, index 1 corresponds to the particle with the highest energy, so that q > 0. The 
longitudinal component, qg, may be either positive, in case the vector difference q lies in the 
pair hemisphere, or negative, in the opposite case. The transverse component, q t , is positive 
definite. 

The experimental three-dimensional correlation functions C are defined, in a small phase 
space volume around each triplet of Qg, Qt sidc and Qt out ( or qt, qt and go) values, as the number 
of like-charge pairs in that volume divided by the number of unlike-charge pairs: 

n _ N w + W + + N n - W - N hkc 

- i *7r+7r _ -'•unlike 

In order to have adequate statistics in each bin, a bin size of 40 MeV was chosen in each 
component of q, which is larger than the estimated detector resolution of 25 MeV [5]. 

Long-range correlations are present in the correlation function C, due to phase space limi- 
tations and charge conservation constraints. In addition, the choice of unlike-sign pairs as the 
reference sample adds further distortions to the correlation function, due to pions from reso- 
nance decays. To reduce these effects, we introduced the (double) ratio C of the correlation 
functions C in the data and in a sample of 7.2 million Jetset 7.4 [23] multihadronic Monte Carlo 
(MC) events, without BECs: 

^DATA atDATA /atDATA 

_ iv likc / iv unlike / a\ 

G " CMC " ^0/^ • W 

The Monte Carlo samples are processed through a full simulation of the OPAL detector [24J. 
The simulation parameters of the generator were tuned in [25] . 

The dependence of the correlation functions C'(q t , qg, go) and C'(Qg, Qt aide , Qt out ) on the 
pair average 4-momentum K has been analyzed by selecting pions in different intervals of two 
components of K: the pair rapidity 



\Y\ = - In 

1 1 2 



(gi+jg) + (j>e,i+Pt,2) 
{E x + E 2 ) - (P£,l + pt,2) 



(5) 



and the pair average transverse momentum with respect to the event thrust direction 



h = \m,i+p % 2)\. (6) 

The differential \Y\ and k t distributions, ■4^ and J^, of the data are shown in Fig. [TJ The 
same distributions for Jetset events are also presented in Fig. [TJ the comparison shows a good 
agreement between data and Monte Carlo events. 
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Figure 1: (a) Histogram of the differential distribution in the pair rapidity \Y\ and (b) in the 
pair mean transverse momentum kt of the data (dots) and Jetset events (line). The number 
of pairs in the Monte Carlo sample has been normalized to the number of pairs in the data 
sample. 

The dependence of C and C on K has been studied in three bins of \Y\ (0.0 < \Y\ < 0.8, 
0.8 < \Y\ < 1.6, 1.6 < \Y\ < 2.4) and five bins of k t (0.1 < k t < 0.2 GeV, 0.2 < k t < 0.3 GeV, 
0.3 < ifet < 0.4 GeV, 0.4 < k t < 0.5 GeV, 0.5 < k t < 0.6 GeV). In this domain, a total of 47.3 
million like-charge and 54.7 million unlike-charge pairs have been analysed. 

3 The experimental correlation functions 

Samples of two-dimensional projections of the correlation function C(Qe, Qt sidc , Qt out ) f° r a 
single bin of \Y\ and k t are shown in Fig. [2] for the data and the MC Jetset events. For 
the example shown @, the bin corresponding to pair rapidities and transverse momenta in the 
intervals 0.8 < \Y\ < 1.6 and 0.3 GeV < k t < 0.4 GeV was chosen. Small (< 0.2 GeV) values 
of Qtout and of Qi nave been required in the (Q t , Qt side ) and in the (Qt Bide , Qt ont ) projections, 
respectively. Bose-Einstein correlation peaks are visible in the data at low Qe,Qt sidc ,Qt out but 
they are not present in the Monte Carlo samples. The same two-dimensional projections for 
the correlation function C'(Qe, Qt sidE , Qt out ) are presented in Fig. [3](a) and (b). Also shown, in 
Fig. Etc), (d) and (e), are the one-dimensional projections for low (< 0.2 GeV) values of the 
other two variables. 

2 Files of the three-dimensional correlation functions will be made available in the Durham HEP database. 
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OPAL 




Figure 2: Two-dimensional projections of the correlation function C(Q£,Q taidc ,Qt ont ) for 0.8 
< \Y\ < 1.6 and 0.3 GeV < k t < 0.4 GeV for the data ((a) and (b)) and for Jetset MC events 
((c) and (d)). Q tout < 0.2 GeV in (a) and (c); Q e < 0.2 GeV in (h) and (d). 

The two-dimensional (qe,qo) arid the one- dimensional q t projections of the correlation func- 
tion C(q t , qt, qo) in the bin 0.8 < \Y\ < 1.6 and 0.3 GeV < k t < 0.4 GeV are shown in Fig. HJ 
for data and Jetset events. Narrow cuts (< 0.2 GeV) on the other variables have been applied 
to make the projections. The combination [(qf + qf) — q^\ of the three variables is an invariant 
greater than zero. This condition and the bound on the pair rapidity constrain the correlation 
function to be different from zero only in a limited region of the (qt, qo) plane, as can be seen in 
Fig. H] (a) and (c). The (qe,qo) and (qe,qt) projections of the correlation function C'(q t ,q£,qo) 
are shown in Fig. [5] together with the one- dimensional projections, for small (< 0.2 GeV) values 
of the other variables. BEC enhancements are clearly seen in both the q t and qi projections, 
Fig. El (c) and (e). Fig. (d), on the other hand, shows that the range available to the variable 
go is quite restricted, and that no Bose-Einstein peak can be observed. 
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Figure 3: Two-dimensional ((a) and (b)) and one-dimensional ((c), (d) and (e)) projections of 
the correlation function C'(Q e ,Q tsidc ,Q tout ) for 0.8 < \Y\ < 1.6 and 0.3 GeV < k t < 0.4 GeV. 
Qt out < 0.2 GeV in (a), Qe < 0.2 GeV in (b). In (c), (d) and (e) the one- dimensional projections 
are obtained for low values (< 0.2 GeV) of the remaining two variables. 



4 Parameterizations of the correlation functions 

To extract the spatial and temporal extensions of the pion source from the experimental corre- 
lation functions, the Bertsch-Pratt (BP) [26] 

C'(Qi, Qt Bidc , Qwt) 

= N(l + Ae~ (Q ^ +Q tdXi dc + ^^ Q ts;dc , Q w ) (7) 

and the Yano-Koonin (YK) [27] 

C'(q u q e , q ) = N(l + Ae-^^ +72fe -^ )2/? ? +72( * ) -^ )2i? o ) )F(g t , q e , q ) (8) 

parameterizations were fitted to the measured correlation functions in all intervals of k t and 
\Y\. 

In both parameterizations, N is a normalization factor while A measures the degree of 
incoherence of the pion sources, and is related to the fraction of pairs that actually interfere. 
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Figure 4: Two-dimensional (qe, qo ) and one-dimensional q t projections of the correlation function 
C(q t , qe, qo) f° r data ((a) and (b)) and Jetset events ((c) and (d)). The correlation function was 
measured in the bin 0.8 < \Y\ < 1.6 and 0.3 GeV < k t < 0.4 GeV. It was required q t < 0.2 
GeV in (a) and (c). In (b) and (d) the one-dimensional projections are obtained for low values 
(< 0.2 GeV) of the remaining two variables. 



The two parameters iV and A, whose product determines the size of the BEC peak, are however 
significantly (anti) correlated: this limits the interpretation of A and the comparison of its values 
between the two parameterizations. 

The two functions F(Q e , Qt sidc , Qt out ) = (l+eion g ^+et sidc <5t sidc +e tout <5t out ) and F(q t , q e , q Q ) = 
(1 + 5 t q t + 5iqe + 5 q ), where e ; and 5\ are free parameters, were introduced in Eq. (7) and 
(8) to take into account residual long-range two-particle correlations, due to energy and charge 
conservation. 

The interpretation of the other free parameters in Eq (7), is the following: 

• -Rt side and -Riong are the transverse and longitudinal source radii in the LCMS, i.e. the 
longitudinal rest frame of the pair; 
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Figure 5: Two-dimensional projections of the correlation function C'(q t , qi, qo): (qi, qo) for q t < 
0.2 GeV in (a) and (qi,q%) for q < 0.2 GeV in (b). One- dimensional projections ((c), (d) and 
(e)) of C'(q t , qe, qo), obtained for low values (< 0.2 GeV) of the remaining two variables. The 
correlation function has been measured in the bin 0.8 < \Y\ < 1.6 and 0.3 GeV < k t < 0.4 
GeV. 



• Rt out and the cross-term -Ri ng,t out are a combination of both the spatial and temporal 
extentions of the source. Under certain assumptions [13], the difference t — R% ., ) is 
proportional to the duration of the particle emission process, and -Rion g ,t out to the source 
velocity with respect to the pair rest frame [22]. 

In the YK function Eq. (8), where 7 = I/a/1 — v 2 , the free parameters are interpreted as 
follows: 

• v is the longitudinal velocity, in units of c, of the source element in the CMS frame; 

• i?o measures the time interval, times c, during which particles are emitted, in the rest 
frame of the emitter (source element). Difficulties in achieving reliable results for the 
time parameter in YK fits have been reported in the literature [2HJ, due to the limited 
phase-space available in 7 2 (go — vqi) 2 ] 
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• R t and Re are the transverse and longitudinal radii, i.e. the regions of homogeneity of 
the source, in the rest frame of the emitter. 

The parameters Rq, R t and Re do not depend on the frame in which the correlation function 
has been measured, since they are evaluated in the rest frame of the source element. 

The two parameterizations are not independent [13], so that a comparison between the BP 
and the YK fits represents an important test. 

5 Results 

Minimum \ 2 fits of the Bertsch-Pratt and the Yano-Koonin parameterizations to the experimen- 
tal correlation functions were performed using the MINUIT [29] program. The error associated 
to each entry of the three-dimensional matrices C and C was computed attributing a Poisso- 
nian uncertainty to the number of like and unlike charge pairs in the corresponding bin. The 
fit range allowed to each variable was set between 40 MeV and 1 GeV. The region below 40 
MeV was excluded to avoid problems of detector resolution and poorly reconstructed or split 
tracks which mimic two like charged particle tracks with very low q. In Sections 5.1 and 5.2 
the results of the fits are presented. Sources of systematic uncertainties on the fit parameters 
are discussed in Section 5.3. Section 5.4 is devoted to a comparison between the BP and the 
YK parameterizations. 

5.1 Bertsch-Pratt fits 

The best-fit parameters of the BP function, Eq. (7), are listed in Tabled], and their dependence 
on \Y\ and k t is shown in Fig. [6] Errors in Fig. [6] include both statistical standard deviations 
as given by the fit program □ and systematic uncertainties (discussed in Section 5.3), added 
in quadrature. One notes that there is only a minor dependence on the rapidity, but some 
parameters depend on k t . In more detail: 

• A varies between 0.25 and 0.4. The coefficient of correlation between the parameters A 
and iV is about —0.35, almost independent of fc t ; 

• -Rt sidc > -^t out an d> l ess markedly, Rf ong decrease with increasing k t . The presence of correla- 
tions between the particle production points and their momenta is an indication that the 
pion source is not static, but rather expands during the particle emission process. R 2 ong 
is larger than the corresponding transverse parameter R% . , in agreement with a pion 
source which is elongated in the direction of the event thrust axis |10j ; 

• the cross-term parameter Rf ong tout is compatible with zero, apart from a few bins at the 
highest rapidity interval. This result may be explained [13] assuming that the source 
velocity, measured with respect to the rest frame of the pion pair, is close to zero; 

• the difference between the "out" and "side" transverse parameters, {R 2 ont — R 2 ) for 
\Y\ < 1.6 is positive at low k t , then it decreases and becomes negative for k t > 0.3 GeV. 
In the highest rapidity interval, 1.6 < |Y| < 2.4, {R 2 — R 2 . ) is compatible with zero, 
for all kt- As a consequence, it is not possible to estimate the particle emission time from 
{R 2 . — R 2 j ); 

\ l<out L sidc ' 1 

3 The HESSE algorithm in MINUIT calculates the error matrix inverting the matrix of the second derivatives 
of the fit function with respect to the fit parameters. 
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Figure 6: Best-fit parameters of the Bertsch-Pratt parameterization, Eq. (7), to the correlation 
function C f (Qe,Qt sidc ,Qt out ), as a function of kt, for different intervals of rapidity \Y\. The 
correlation functions were measured in the LCMS frame. Horizontal bars represent bin widths 
and vertical bars include both statistical and systematic errors, (a ) the normalization factor N; 
(b) the incoherence parameter X; (c) the cross term R\ ong tout ; (d) the parameter R^ out ; (e) the 
squared longitudinal correlation length Rf ong and (f) the squared transverse correlation length 

L sidc 

• the parameters e Y are not negligible: the function F(Qg, Qt side , Qt out ) typically differs from 
unity for at most 15% to 20% at Q- x w 1 GeV. 



5.2 Yano-Koonin fits 

Tableland Fig. [7] show the parameters of the YK fits, Eq. (8), in different \Y\ and k t intervals. 
Error bars in Fig. [7| include both statistical and systematic uncertainties, added in quadrature. 
It can be seen that: 
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Figure 7: Best-fit parameters of the Yano-Koonin parameterization, Eq. (8), to the correlation 
function C'(q t , qe, qo), as a function of k t , for different intervals of rapidity \Y\. The correlation 
functions were measured in the event centre-of-mass frame. Horizontal bars represent bin widths 
and vertical bars include both statistical and systematic errors, (a) the normalization factor 
N; (b) the parameter A; (c) the source velocity v; (d) the time parameter R^; (e) the squared 
longitudinal correlation length Rj and (f) the squared transverse correlation length R^. 



• the parameter A is almost independent of rapidity and increases with k t , reaching values 
of about 0.5 for the largest k t values. It is however significantly anticorrelated with the 
parameter N, the correlation coefficient increasing in absolute value from about —0.50 at 
low kt up to -0.80 for k t > 0.4 GeV; 

• both R% and R\ decrease with increasing k t and The longitudinal radii are larger 
than the transverse radii. This agrees with an expanding, longitudinally elongated source; 

• Rq is compatible with zero at high rapidities, and assumes negative values for |Y| < 1.6. 
This excludes an interpretation of Rq/c in terms of the time duration of the particle 
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emission process; 



• those of the parameters 8[ which are not negligible, contribute typically 10% to 15% to 
the function F(q t ,q iy q ) at large g ; ; 

• the source velocity v does not depend on k t , but it is strongly correlated with the pair 



The dependence of v on \Y\ can also be presented [301 H] in terms of a plot a la GIBS, i.e. the 
Yano-Koonin rapidity 



as a function of the pair rapidity |Y|. Yyk measures the rapidity of the source element with 
respect to the centre-of-mass frame: a non-expanding source would therefore correspond to 
Yyk ~ for any |Y|. On the other hand, for a boost-invariant source 0, the strict correlation 
Yyk = \y\ is expected [T3| |2~T|. since only the source elements which move with velocities close 
to the velocity of the observed particle pair contribute to the correlation function. In Fig. [8] 
the Yano-Koonin rapidity Yyk is shown as a function of the pair rapidity. Since in a given 
\Y\ interval the parameter v is almost independent of k t , see Fig. [7(c), each Yyk is computed, 
according to Eq. (9), using the average value of v over all k t in that \Y\ bin. Each |Y| has been 
computed as the weighted average of the corresponding bin, rather than the centre of the bin. 
A clear positive correlation between Yyk and |Y| is observed, even if Yyk < |Y| at the largest 
pair rapidities. This is in agreement with a pion source which is emitting particles in a nearly 
boost-invariant way. 

To try to understand the YK fit results of the parameter it is useful to analyse the two- 
dimensional projection (qe,qo) of the correlation function C'(q t ,qi,qo) after the longitudinal 
boost to the rest frame of the source element. We then introduce g^ oost = ^(^ — vqo) and 
gboost _ _ V q £ ^ where the best-fit parameter v is used to boost the variables. In Fig. 0(a) 
the two-dimensional (|^ oost |,|go° ost |) projection of C is presented. The phase space available to 
jgboostj j g ij m ited, when q t approaches 0, and the one- dimensional |gQ OOSt | projection (Fig. El(b)) 
is approximately flat: it is not possible to distinguish any peak due to Bose-Einstein correlations 
and, for most rapidity and k t intervals, the fitted Rq have negative values. In Fig. [9](b), the 
solid line shows the one-dimensional |gQ OOSt | projection (|g^ oost | < 0.2 GeV, q t < 0.2 GeV) of 
the YK fit, Eq. (8); the line is an increasing function of gQ OOSt , because of the negative value of 
-Rq. Similar limitations in the temporal acceptance have been reported in the literature [28J. 
On the other hand, the \q^ oost \ projection for |g£ oost | < 0.2 GeV and q t < 0.2 GeV, Fig. [DTc), 
shows a clear BEC peak at small |g boost |, reproduced by the one dimensional |g£ OOSt | projection 
of the best-fit YK function (solid line). 

5.3 Systematic effects 

The systematic uncertainties of the fit parameters and the stability of the results concerning 
the dependence of the transverse and longitudinal radii on k t was studied by considering a 
number of changes with respect to the reference analysis. The following changes were taken 
into account: 

4 A source expands boost-invariantly in the longitudinal direction if the velocity of each element is given by 
v = z/t, where t and z are, respectively, the time elapsed since the collision and the longitudinal coordinate 
of the element, in the centre-of-mass frame. In that case, particle emission happens at constant proper times 



rapidity. 




(9) 
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Figure 8: The Yano-Koonin rapidity Iyk plotted versus the pion pair rapidity \Y\. Each \Y\ 
was computed as the weighted average of the corresponding bin. Yyk values were computed by 
means of Eq. (9), using the average value of v over all k t in that \Y\ bin. Horizontal bars are 
r.m.s. deviations from the average. Vertical bars include both statistical and systematic errors. 
Also shown is the line Yyk=\Y\, corresponding to a source which expands boost-invariantly. 

• A correction was applied to the correlation functions, based on the Gamow factors [31J, 
in order to take into account final-state Coulomb interactions between charged pions. 

• The analysis was repeated with more stringent cuts in the selection: a maximum momen- 
tum of 30 GeV instead of 40 GeV and a charge unbalance smaller than 0.25 per event 
instead of 0.4. 

• The fits were repeated changing the upper bound of the fit range from 1 GeV to 0.8 GeV. 

In the cases listed above, we found negligible differences in the parameters with respect to the 
reference analysis. The systematic effect on the correlation function C, due to the Monte Carlo 
modelling, was assumed negligible. 

• The correlation functions were measured in bins of 60 MeV, instead of 40 MeV, to test the 
stability of the fits. Bin widths larger than 60 MeV would prevent a correct reconstruction 
of the BEC peak, which is about 300-^400 MeV wide. 

• Possible non-Gaussian shapes of the correlation functions at low q were tested replacing 
the Gaussian functions in the BP and YK parameterizations with first order Edgeworth 
expansions [32] of the Gaussian. The x 2 /DoF of the two fits were found to be comparable. 
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Figure 9: (a) The two-dimensional projection (\qi° ost \, \Qo° ost \), after the longitudinal boost to 
the source element rest frame, measured for pion pairs in the rapidity interval 0.8 < \Y\ < 1.6 
and with mean transverse momenta in the range 0.3 GeV < k t < 0.4 GeV. The projection was 
made requiring q t < 0.2 GeV. 

(b) The one- dimensional projection in |gQ OOSt | (]g]? oost | < 0.2 GeV). The curve is the one- 
dimensional projection of the Yano-Koonin three-dimensional best-fit function. 

(c) The one-dimensional projection in |gj? oost | (|qo oost | < 0.2 GeV). The curve is the one- 
dimensional projection of the Yano-Koonin three-dimensional best-fit function. 



Systematic errors on the fit parameters have been computed adding in quadrature the deviations 
from the standard fit; they are reported in Tables [1] and [2j 

Assuming simple linear dependences of the squared BP and YK longitudinal and transverse 
radii on k t , we measured the slopes, dRf/dk t , by minimum x 2 fits. Fits were performed on 
the radii of the reference analysis, with statistical errors only. The systematic errors on the 
slopes were then estimated comparing the slopes from the reference analysis with the slopes 
from the systematic checks listed above. Table [3] shows the best-fit slopes with errors. In all 
cases a decrease of the radii with increasing k t is favoured even if, in one rapidity interval, the 
longitudinal BP radius is compatible with independence on k t . 

To investigate further the decrease of the radii on k t , the YK and BP functions were fitted 
to the correlation function C, Eq. (1). Larger (about 30%) squared transverse and longitudinal 
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radii with respect to the correlation function C are obtained in this case. However, the slopes of 
the linear dependences of the squared radii on k t are the same, within uncertainties, for C and 
C . A comparison of the YK best-fit parameters from minimizing \ 2 values and from maximizing 
a likelihood function [33] has been done for the correlation function C. The differences between 
the parameters fitted with the two techniques were negligible. 

One more check was done on the YK transverse radius R t : we computed the one-dimensional 
projection C'(q t , 0, 0) of the three-dimensional correlation function C'(q t , qi, q ), by requiring qi 
and q < 0.08 GeV, and we fitted the function 

C'{q t ) = N(l + \e- q ? R <) (10) 

to the projection. We first checked that the best-fit R\ is compatible, within errors, to the one 
we obtain if the right-hand side of Eq. (10) is multiplied by a "long-range" factor (1 + 5 t <Zt)- 
Based on the same one-dimensional projection C'(qt,0, 0), we also measured the transverse 
correlation length in a fit-independent way [31] , introducing the parameter R t 



6 _ ' wW /„2, Jft 2 [C'fa,0,0)-l]d9. , , 



i.e. the inverse variance of the correlation function for small qt values 0. We found that R t , 
computed using Eq. (11), agrees with the best-fit R t from Eq. (10); the slope of the linear 
decrease is about 20% smaller than the one measured with three-dimensional YK fits, Eq. (8). 

The standard analysis was also repeated for a subsample of events classified as two-jets by 
the Durham jet-finding algorithm [35]. The subsample was defined by setting the resolution 
parameter at y cut = 0.04. The dependences of the best-fit parameters on \Y\ and k t are similar to 
those found for the inclusive sample of events. In particular, the longitudinal and the transverse 
radii decrease with increasing k t . However, the radii measured in the case of two-jet events are 
smaller, by about 10%, than in the inclusive sample [5]. An increase of the "jettyness" of the 
two-jet subsample, obtained using a smaller y cnt (y cut = 0.02) in the jet-finding algorithm, does 
not change significantly the behaviour of the parameters. 

5.4 Comparison between BP and YK fits 

The following relations should hold between the correlation lengths of the BP and YK functions 
measured in the LCMS and CMS frames, respectively 



R h = R i 2 2 (12) 

-^long = 7LCMs(-^£ + AxMS-^o) (13) 

("^tout _ R t M J = A 2 7lcms(-Ro +PtcMS R i)- ( 14 ) 
In Eq. (13) and (14) /3lcms is the velocity of the source element measured in the LCMS, i.e. with 
respect to the pair longitudinal rest frame; 7lcms = — /^lcms- ^ n Eq- (14) fi\ = ^ E f+ E2 ) , 

where the brackets stand for the average over all pion pairs in the given \Y\ and k t range. For 
a boost-invariant source, Plcms — and Eqs. (13) and (14) reduce to: 

Aong = R £ (15) 

(RL.-KJ = AX (i 6 ) 



5 In the actual estimate of (q?) we have computed ^ - vl^, ( - gt ' 0,0 ^ — J. where N is the normalization parameter 

X *' 2j[C'(g t ,0,0)-Af] ' 1 

of the fit Eq. (10) and each q t has been taken as the central value of the corresponding 40 MeV bin. 
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Figure 10: (a) (d) (g) The best-fit longitudinal radius Rf ong of the Bertsch-Pratt parameteriza- 
tion (open dots) compared with the Yano-Koonin longitudinal radius R\ (full dots), (b) (e) (h) 
The BP transverse correlation length R% . (open dots ) compared with the YK transverse cor- 
relation length R% (full dots), (c) (f) (i) The difference of the BP transverse radii (R% —Rt ) 
(open dots) compared with the YK time parameter R^ times 01 (full dots). Errors on the 
parameters include both statistical and systematic uncertainties, added in quadrature. 

In Fig. [10] the best-fit BP parameters -R^g, Rt sidc an d (Rt out ~ ^t id ) are compared with the 
YK parameters Rj, R% and /9 t 2 Ro- 

The longitudinal parameter -R^g is systematically larger than Rj in all the rapidity intervals 
analyzed (Fig. [TOT a). (d) and (g)). According to Eq. (13), -Rj^g > Rj corresponds to /?lcms 
greater than zero, in agreement with a pion source whose expansion is not exactly boost- 
invariant. 

The equality of the transverse parameters -Rt sidc an d R% > Eq. (12), is confirmed within errors, 
with possible deviations at low k t (Fig. fTQTb). (e) and (h)). 

The negative values of Rq and (R% out — R% . d ) appearing in the two first rapidity intervals 
(Fig. [TOTc). (f) and (i)) prevent an interpretation in terms of the time duration of the particle 
emission process. Negative values of Rq have been suggested |36j as possible indicators for 
opacity of the source, i.e. surface dominated emission. A dependence of (R% ont — R% id ) on 
k t similar to the one shown in Fig. [TOTc) and (f) has been reported in heavy- ion collision 
experiments [37] . 
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6 Conclusions 



An analysis of Bose-Einstein correlations in e + e~ annihilation events at the Z° peak performed 
in bins of the average 4-momentum of the pair, K, has been presented for the first time. Based 
on this, dynamic features of the pion emitting source were investigated. Previous BEC analyses, 
not differential in K, were not sensitive to these features. 

Using the Yano-Koonin and the Bertsch-Pratt formalisms, the correlation functions were 
studied in intervals of two components of K: the pion pair rapidity \Y\ and the mean transverse 
momentum k t . We found that the transverse and longitudinal radii of the pion sources decrease 
for increasing k t , indicating the presence of correlations between the particle production points 
and their momenta. The Yano-Koonin rapidity scales approximately with the pair rapidity, 
in agreement with a nearly boost-invariant expansion of the source of pions. Limitations in 
the available phase space did not allow measurement of the duration of the particle emission 
process. 

Similar results have been observed in more complex systems, such as the pion sources cre- 
ated in pp and heavy-ion collisions, which are now complemented with such measurements in 
the simpler hadronic system formed in e + e~ annihilations. 
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0.0 < \Y\ < 0.8 



0.1 < fc t < 0.2 GeV 



0.2 < fc t < 0.3 GeV 



0.3 < fc t < 0.4 GeV 



0.4 < fc t < 0.5 GeV 



0.5 < fc t < 0.6 GeV 



N 
A 

<» t (fm 2 ) 

*L. ( fm2 ) 

*Lg,w (fm 2 ) 
et out (GeV- 1 ) 
et Bldc (GeV- 1 ) 
eiong (GeV' 1 ) 



0.974 ±0.003 ±0.057 
0.286 ±0.011 ±0.067 
0.60 ±0.07 ±0.18 
0.50 ±0.03 ±0.14 
1.09 ±0.11 ±0.37 
-0.06 ±0.08 ±0.14 
-0.091 ±0.004 ±0.060 
-0.123 ±0.004 ±0.071 
0.081 ±0.004 ±0.023 



0.996 ±0.004 ±0.042 
0.364 ±0.009 ±0.061 
0.36 ±0.03 ±0.10 
0.38 ±0.02 ±0.11 
0.72 ±0.04 ±0.17 
0.020 ±0.036 ±0.037 
-0.056 ±0.004 ±0.035 
-0.130 ±0.004 ±0.061 
0.048 ±0.004 ±0.017 



1.011 ±0.004 ±0.040 
0.429 ±0.012 ±0.047 
0.294 ±0.020 ±0.079 
0.37±0.02±0.11 
0.82 ±0.05 ±0.16 
-0.065 ± 0.028 ±0.031 
-0.037 ±0.004 ±0.027 
-0.140 ±0.004 ±0.071 
0.019 ±0.005 ±0.015 



1.003 ±0.007 ±0.040 
0.398 ±0.013 ±0.044 
0.174 ±0.013 ±0.050 
0.30 ±0.02 ±0.10 
0.60 ±0.04 ±0.15 
-0.023 ± 0.022 ±0.007 
-0.016 ± 0.005 ±0.018 
-0.18 ±0.01 ±0.10 
0.016 ±0.007 ±0.019 



1.016 ±0.009 ±0.052 
0.337 ±0.016 ±0.063 
0.169 ±0.014 ±0.051 
0.22 ±0.03 ±0.10 
0.75 ±0.06 ±0.18 
-0.121 ± 0.024 ±0.065 
-0.003 ±0.006 ±0.018 

-0.24 ±0.01 ±0.13 
-0.018 ±0.008 ±0.038 



x 2 /DoF 



16389/15617 



16080/15617 



15596/15617 



15864/15617 



15439/15617 



0.8 < \Y\ < 1.6 



N 
A 

RL t (fm 2 ) 
(fm 2 ) 

Rln S (fm 2 ) 
*U,W (fm 2 ) 
ew (GeV- 1 ) 
et sldo (GeV- 1 ) 
eiong (GeV' 1 ) 



0.972 ± 0.003 ± 0.049 
0.315 ±0.008 ±0.070 
0.62 ±0.04 ±0.20 
0.52 ±0.03 ±0.13 
1.06 ±0.08 ±0.35 
0.019 ± 0.055 ±0.076 
-0.070 ±0.004 ±0.049 
-0.106 ±0.004 ±0.059 
0.066 ±0.004 ±0.026 



0.990 ±0.004 ±0.075 
0.386 ±0.008 ±0.064 
0.38 ±0.02 ±0.11 
0.38 ±0.02 ±0.11 
0.71 ±0.04 ±0.16 
-0.029 ±0.026 ±0.036 
-0.046 ±0.004 ±0.035 
-0.104 ±0.004 ±0.051 
0.035 ±0.004 ±0.026 



1.017 ±0.005 ±0.052 
0.393 ±0.011 ±0.053 
0.271 ±0.016 ±0.079 
0.34 ±0.02 ±0.10 
0.65 ±0.04 ±0.15 
-0.036 ±0.023 ±0.025 
-0.033 ±0.004 ±0.033 
-0.131 ±0.005 ±0.064 
-0.003 ±0.005 ±0.036 



1.019 ±0.007 ±0.057 
0.379 ±0.013 ±0.055 
0.204 ±0.014 ±0.062 
0.272 ±0.021 ±0.081 

0.64 ±0.05 ±0.16 
-0.061 ± 0.022 ±0.045 
-0.015 ± 0.005 ±0.027 
-0.161 ± 0.005 ±0.088 
-0.028 ±0.006 ±0.047 



1.024 ±0.010 ±0.066 
0.318 ±0.016 ±0.062 
0.141 ±0.015 ±0.053 
0.226 ±0.026 ±0.079 

0.50 ±0.05 ±0.12 
-0.034 ±0.021 ±0.042 
0.007 ±0.007 ±0.009 
-0.23 ±0.01 ±0.14 
-0.060 ±0.009 ±0.072 



x/DoF 



15856/15617 



15745/15617 



15658/15617 



15895/15617 



15592/15617 



1.6 < \Y\ < 2.4 



N 
A 

<» t (fm 2 ) 
(fm 2 ) 
(fm 2 ) 
(fm 2 ) 
ew (GeV- 1 ) 
et sldo (GeV- 1 ) 
eiong (GeV- 1 ) 



^t 2 - d 

L sid 

p2 
^Mon. 



0.991 ±0.003 ±0.082 
0.261 ±0.008 ±0.079 
0.54 ±0.04 ±0.19 
0.53 ±0.03 ±0.15 
1.13 ±0.09 ±0.41 
-0.05 ±0.05 ±0.13 
-0.102 ±0.004 ±0.070 
-0.134 ±0.004 ±0.079 
0.045 ±0.004 ±0.056 



1.019 ±0.005 ±0.069 
0.307 ±0.008 ±0.072 
0.35 ±0.02 ±0.11 
0.34 ±0.02 ±0.10 
0.61 ±0.04 ±0.17 
0.012 ±0.029 ±0.033 
-0.063 ±0.004 ±0.048 
-0.130 ±0.004 ±0.068 
-0.003 ±0.005 ±0.045 



1.066 ±0.005 ±0.078 
0.299 ±0.011 ±0.065 

0.35 ±0.03 ±0.10 
0.279 ±0.023 ±0.077 

0.78 ±0.06 ±0.19 
-0.137 ± 0.033 ±0.076 
-0.060 ±0.005 ±0.043 
-0.167 ±0.005 ±0.082 
-0.046 ±0.005 ±0.053 



1.055 ±0.008 ±0.074 
0.264 ±0.014 ±0.074 
0.219 ±0.017 ±0.064 
0.229 ±0.026 ±0.072 

0.62 ±0.05 ±0.17 
-0.148 ± 0.024 ±0.077 
-0.027 ±0.006 ±0.028 

-0.19 ±0.01 ±0.11 
-0.078 ±0.007 ±0.072 



1.07 ±0.01 ±0.10 
0.24 ±0.02 ±0.11 
0.25 ±0.03 ±0.11 
0.169 ±0.034 ±0.085 
0.54 ±0.07 ±0.13 
-0.09 ±0.04 ±0.11 
0.16 ±0.01 ±0.12 
-0.31 ±0.01 ±0.17 
-0.15 ±0.01 ±0.13 



x7DoF 



15966/15617 



15866/15617 



15735/15617 



15235/15617 



15279/15617 



Table 1: Results of the Bertsch-Pratt fits, Eq. (7), to the experimental three-dimensional correlation functions C'(Q£,Q ts . dc ,Q tout ) over 
the range 0.04 < Qg, Qt Bide , Qt out < 1-0 GeV. The first errors are statistical and the second systematic. The quality of the fits is indicated 
by the value of x 2 /DoF, which ranges from 0.98 to 1.05. 



0.0 < \Y\ < 0.8 


0.1 < fct < 0.2 CjreV 


0.2 < fc t < 0.3 Lrev 


0.3 < fc t < 0.4 LreV 


0.4 < fct < 0.5 CjreV 


0.5 < fct < 0.6 CjceV 


N 




A AA o I A AAO 1 A A1 A 

0.993 ± 0.003 ± 0.010 


1 A A 7 I A AAO I A AAA 

1.004 ± 0.003 ± 0.008 


AAAP" 1 A A A 7 1 A A 7 ^7 

0.985 ± 0.004 ± 0.017 


A A 7 A 1 A A A f 1 A A 7 7 

0.946 ± 0.005 ± 0.041 


A Ci f 1 A A 7 1 Aid 

0.85 ± 0.01 ± 0.14 


A 




A A A A 1 A A1 A I A A/lf7 

0.266 ± 0.012 ± 0.067 


A o h» A | A AAA I A A A 

0.358 ± 0.009 ± 0.056 


A A A 7 1 A A 7 A 1 A AO ^ 

0.441 ± 0.012 ± 0.035 


A A A A 1 A A 7 O 1 A AAO 

0.440 ± 0.013 ± 0.023 


A 7 O A 1 A A 7 ^7 1 AAAA 

0.482 ± 0.017 ± 0.069 






a 7 a i a 1 a i a o n 

0.10 ± 0.19 ± 0.32 


a nnn i a a 7 a i a ai n 

0.288 ± 0.042 ± 0.018 


A OAA 1 A AOA 1 A A7 O 

0.320 zt 0.030 zt 0.013 


A A 7 A 1 A A O 7 1 A A 7 A 

0.249 ± 0.034 ± 0.010 


A A 7 1 I A A07 I A A07 

0.2H ± 0.031 ± 0.031 


Rl (fm 2 ) 




n ^0 -U n on -i_ n 1 ^ 

— U.OZ zt U.ZU it U.10 


— U. 154 zt U.U44 zt U.U / 


n ooa -U n noo -u n hqq 
— U.ZzU zt U.Uzz zt U.Uoo 


0A7 4- A AI 7 4- A A81 


A 1 1 A 4- A AAO 4- A A^O 

— u.iiu ± u.uuy ± u.uoo 


Rl (fm 2 ) 




7^ + 04 + 1 7 

yj. 1 u —I— \J ,\J^± _l_ U.l 1 


£,0 + 09 + 019 

yj.'jyj —1— u.UZi i yj . A. Li 


n 41 _l n no + o 1 


717 4-0 011 + 081 


1 07 4- 007 + 048 

U.lfO —1— yJ .UU f —1— U.U4:G 


Rj (fm 2 ) 




1.15 ±0.15 ±0.34 


0.74 ± 0.04 ± 0.13 


0.66 ± 0.03 ± 0.11 


0.46 ± 0.02 ± 0.10 


0.31 ±0.02 ±0.11 


So (GeV- 1 ) 
it (GeV" 1 ) 




-0 045 ± 007 ± 069 

V 7 . V / 1 J — 1 — V 7 . V 7 V 7 f — 1 — V / ■ V 7 V 7 - 7 


-0 014 ± 006 ± 052 

V / ■ V 7 1 1 — 1 — V / ■ V 7 V 7 V 7 — 1 — V / ■ \ > 7 


008 ± 007 ± 019 

v 7 . v / v / v. 7 — 1 — v / ■ v 7 v 7 j — 1 — v 7 . v / 1 - / 


061 ± 008 ± 062 

v / ■ v 7 v 7 1 — 1 — v / ■ v 7 v 7 l ' — 1 — v 7 . v / \ 7 —- 


156 ± 011 ± 091 

v / ■ 1 • 7 v 7 — 1 — v 7 . v / 1 1 — 1 — v 7 ■ v / - / 1 




-0 099 ± 004 ± 029 

V 7 . V / - 7 - 7 — 1 — V 7 • V 7 V 7 1 — 1 — V / ■ V 7 • 7 


-0 089 ± 005 ± 023 

u . uu kj — i — v / ■ \ i \ > ■ j — i — \ i • \ i ^ * j 


-0 065 ± 006 ± 030 

v / ■ v 7 v / . 7 — 1 — v / ■ v 7 v 7 v / — 1 — v/.v/t7V/ 


-0 066 ± 008 ± 036 

v 7 . v / v 7 v 7 — 1 — v 7 ■ v / v / v.7 — 1 — v 7 . v / 1 7 v 7 


-0 054 ± 013 ± 019 

U . U U 1 — 1 — V 7 ■ V / 1 t 7 — 1 — V / ■ V 7 1 -7 


Si (GeV- 1 ) 




038 ± 002 ± 098 


014 ± 002 ± 090 


005 ± 003 ± 051 


—0 001 ± 003 ± 025 


—0 002 ± 005 ±0012 

V 7 . V / V / — w 1 U.UUU 1 U . U -L £1 


x 2 /DoF 


1 3583 /1 1 677 


16008/14375 


1 7555 /1 6338 


1 81 66/1 7554 

IOIUU / J- / OO^l 


1 8709 /I 8939 

IO f yj Li j IO LiKJ Li 


0.8 < \Y\ < 


1.6 












N 




a r\ ,i a i a AAr i a a a 7 

0.948 ± 0.005 ± 0.021 


A t~\ >—?>—! 1 AAA/"' 1 A A A /"* 

0.977 ± 0.006 ± 0.006 


A A r> A 1 A AAA 1 A A 7 A 

0.964 ± 0.009 ± 0.010 


A A 7 7 I A A 7 A I A A P"* 7 

0.914 ± 0.012 ± 0.051 


A A r* A 1 A AAA 1 A ATA 

0.859 ± 0.022 ± 0.072 


A 




a oa \ i a aaa i a a a r* 

0.324 ± 0.008 ± 0.065 


AOAA 1 A A 7 A 1 AAP'^T 

0.380 ± 0.010 ± 0.057 


A 7 A f 1 A A 7 O 1 A AO f 

0.425 ± 0.013 ± 0.035 


A A A A 1 A A 7 A 1 A AAO 

0.464 ± 0.019 ± 0.023 


A A A A 1 A A O A 1 A AAA 

0.464 ± 0.034 ± 0.028 






0.754 ± 0.022 ± 0.061 


a rrAA 1 A A 1 A I A A O A 

0.782 ± 0.014 ± 0.036 


A >~? A A 1 A AI 7 1 A A A !~? 

0.742 ± 0.017 ± 0.027 


A ryryrj \ A A 1 A I A A0 1 

0.777 ± 0.014 ± 0.031 


A f7 A 1 A AAO 1 A A /| A 

0.743 ± 0.023 ± 0.040 


\l (fm 2 ) 




n 1 87 4- A A^7 4- A A7fi 
— U.lof ±U.U04±U.U/0 


n iH/i -u n no a _i_ n n^i 

— U.1U4 zt U.UZ4 zt U.UOl 


n 11/1 4-D nifi-u n n^/i 

— U.114 zt U.U1D zt U.U04 


n 1 no j- n m q -U n n^i 
— U.lUz zt U.Ulo zt U.UOl 


n 1 r\fi -U n m 1 -u n n^o 
— U.1UU zt U.Ull zt u.uoz 


i? t 2 (fm 2 ) 




0c;«4-0 09-l-01 c i 


+ 01 +011 

U.Of ZIZ U.U± ZIZ \J.±± 


n 00 _l n ni _l n 1 n 
xj.OLi zn u.ui zn u. iu 


n 00c _l ni n _l n (-171 
\j.ZiOO zn u.uiu zn u.ui ± 


n 1 c\A + 00Q + 04Q 
U.1U4: zn u.uUiv zn u.u^tiy 


Rj (fm 2 ) 




0.83 ± 0.06 ± 0.27 


0.58 ± 0.03 ± 0.13 


0.46 ± 0.03 ± 0.10 


0.43 ± 0.03 ± 0.12 


0.294 ± 0.032 ± 0.085 


(Jo (GeV- 1 ) 
<$t (GeV" 1 ) 




-0 07 ± 01 ± 10 

U . U 1 — 1— U ■ U 1 _l_ U . J- U 


00 ± 01 ± 12 

V 7 . V / V 7 1 V / ■ V 7 1 1 \_f ■ _L t-t 


071 ± 010 ± 053 


106 ± 012 ± 076 


125 ± 019 ± 084 

U. luU 1 V 7 . V / 1 • J 1 V 7 . V / V 7 1 




—0 068 ± 007 ± 041 


—0 075 ± 008 ± 033 


—0 069 ± 011 ± 043 


—0 046 ± 01 5 ± 036 

\J ,\_)i:\J _1_ U.UIU —1— U.UUU 


—0 079 ± 027 ± 079 

U.UIU —1— U.Uti 1 —1— U.UI tj 


Si (GeV- 1 ) 




099 ± 010 ± 078 

V / ■ V ' ' / ' 7 — 1 — V 7 . V / 1 V / — 1 — V / ■ V 7 J l_ 7 


022 ± 008 ± 099 

V7 ■ V -j -j — 1 — V / ■ V 7 V 7 L ' — 1 — V / ■ V 7 - 7 - 7 


-0 032 ± 008 ± 062 

V / ■ V 7 • 7 — 1 — V / ■ V ' V 7 1. ' — 1 — V 7 . V / \ 7 


-0 058 ± 010 ± 047 

U.UUU — 1 — V 7 ■ V / 1 \ ' — 1 — u . u 1 1 


-0 077 ± 014 ± 060 

U.UI 1 — 1 — V 7 ■ V / 1 1 — 1 — U.UUU 


X 2 /DoF 


8624/71 39 


9778/8788 

t/( 10/0 100 


1 1004/9870 


1 1 365/1051 8 


1 1603/10885 

-LJ-UUu / -LuOOu 


1.6 < \Y\ <■ 


2.4 












N 




/"I n A A I A A7 /"i 1 A Af < 

0.899 ± 0.016 ± 0.054 


A A f O 1 A A A A 1 A A /"* < 

0.963 ± 0.020 ± 0.064 


AAAA 1 AAA1 I A A A A 

0.902 ± 0.021 ± 0.040 


AAAA I AAAA 1 A A P" A 

0.888 ± 0.028 ± 0.052 


A 7 Cl I A A 7 I A OA 

0.48 ± 0.01 ± 0.39 


A 




A O 7 A 1 A A1 A | A AT A 

0.342 ± 0.019 ± 0.078 


A Oh" A 1 A AAA 1 A A^TA 

0.354 ± 0.022 ± 0.072 


A A f A I A A 7 A I A A 7 7 

0.454 ± 0.012 ± 0.041 


A 70A 1 A AOA 1 A AOA 

0.438 ± 0.038 ± 0.030 


1 A /"1 I A A 7 I AAO 

1.26 ± 0.04 ± 0.63 






A n A O I A A1 A I A A A A 

0.893 ± 0.012 ± 0.044 


A A O 7 1 AAAA 1 A A 7 >~7 

0.931 ± 0.008 ± 0.047 


A A AW 1 A AAA 1 A AAO 

0.927 ± 0.009 ± 0.033 


A A 7 A 1 A A 7 A 1 A A 7 A 

0.912 ± 0.012 ± 0.042 


A A O I A A 7 I A 1 1 

0.93 ± 0.04 ± O.ll 


Rl (fm 2 ) 




A AQ1 _l_ A A/1Q -1- A non 

U.Uol ± U.U40 ± u.uzy 


A AAft 4- A HQO 4- A AOO 

U.UUO ± U.Uoz ± v.uzz 


A AI K 4_ A HOfi 4- A nor; 
U.UlO ± U.UZD ± U.Uzo 


A AOA 4- A AQA 4- A AOA 
U.UZU ± U.UoU ± U.Uzo 


A AQ7 4- A AO/1 4- A A1 O 

— U.U04 ± U.UZ4 ± u.uiy 


J? t 2 (fm 2 ) 




70 4-0 09 + 019 
u.oy in u.uz in u.iz 


9Q1 4-00174-0 087 


1 79 4- 01 1 4- 0^,9 


1 KQ 4- 01 A 4- 048 


071 4- 00^, 4- 091 
U.Uf 1 41 U.UUO 41 U.UZ1 


i?, 2 (fm 2 ) 




0.42 ±0.05 ±0.16 


0.42 ±0.05 ±0.15 


0.34 ±0.03 ±0.10 


0.38 ±0.06 ±0.10 


0.12 ±0.03 ±0.11 


<5o (GeV- 1 ) 
St (GeV- 1 ) 




0.03 ±0.04 ±0.12 


0.060 ±0.029 ±0.080 


-0.014 ±0.033 ±0.021 


0.196 ±0.031 ±0.075 


0.00 ±0.03 ±0.16 




-0.033 ±0.023 ±0.082 


-0.087 ± 0.024 ±0.093 


-0.030 ±0.009 ±0.050 


-0.070 ± 0.032 ±0.055 


0.216 ±0.031 ±0.071 


St (GeV- 1 ) 




-0.001 ±0.036 ±0.066 


-0.049 ±0.028 ±0.081 


0.025 ±0.030 ±0.054 


-0.162 ± 0.029 ±0.065 


0.02 ±0.03 ±0.16 


X 2 /DoF 


4168/3804 


5110/4648 


5952/5169 


5775/5490 


5876/5642 



Table 2: Results of the Yano-Koonin fits, Eq. (8), to the experimental three-dimensional correlation functions C'(qt,qe,qo) over the 
range 0.04 < qt,qe,qo < 1-0 GeV. The hrst errors are statistical and the second systematic. The quality of the fits is indicated by the 
value of x 2 /DoF, which ranges from 1.03 to 1.20. 





Jr>± ] 

dRf ong /dk t 
(fnr/GeV) 


radii 

dRlJdh 

(fnr/GeV) 


X XV 

dR\ / dk t 
(fnr/GeV) 


radii 

dR 2 Jdk t 
(far/GeV) 


\Y\ 


< 0.8 


-0.46 ±0.20 ±0.35 


-0.59 ±0.08 ±0.19 


-1.60 ±0.13 ±0.38 


-1.14 ±0.05 ±0.23 


0.8 < 


\Y\ < 1.6 


-0.91 ±0.18 ±0.30 


-0.66 ±0.08 ±0.15 


-1.04 ±0.12 ±0.23 


-0.84 ±0.04 ±0.15 


1.6 < 


|y| < 2.4 


-0.64 ±0.21 ±0.36 


-0.80 ±0.09 ±0.28 


-0.82 ±0.13 ±0.17 


-0.70 ±0.04 ±0.20 



Table 3: Slopes of the linear fits to the dependence of the longitudinal and transverse squared radii of the BP and YK parameterizations 
on k t . Input to the fits are the measured values of Rf ong , -Rt sidc , R\ and R\, reported in Tables^ and\2i The first errors are statistical 
and the second systematic. 



